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ABSTRACT 



We report the discovery of a young planetary-mass brown dwarf in the p 
Oph cloud core. The object was identified as such with the aid of a 1.5-2.4 
/iin low-resolution spectrum obtained using the NIRC instrument on the Keck I 
telescope. Based on the COND model, the observed spectrum is consistent with a 
reddened (A v ~ 15 — 16) brown dwarf whose effective temperature is in the range 
1200-1800 K. For an assumed age of 1 Myr, comparison with isochrones further 
constrains the temperature to ~ 1400 K and suggests a mass of ~ 2-3 Jupiter 
masses. The inferred temperature is suggestive of an early T spectral type, which 
is supported by spectral morphology consistent with weak methane absorption. 
Based on its inferred distance (~ 100 pc) and the presence of overlying visual 
absorption, it is very likely to be a p Oph cluster member. In addition, given 
the estimated spectral type, it may be the youngest and least massive T dwarf 
found so far. Its existence suggests that the initial mass function for the p Oph 
star-forming region extends well into the planetary-mass regime. 

Subject headings: stars: low-mass, brown dwarfs — stars: pre-main sequence - 
infrared: stars 
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Introduction 



Free -floating planetary-mass obj ects have been purpo rted in young clus ters such as a 



Orionis (jZapatero Osorio et al. 



2006| ) . These are 



20021 ) and the Trapezium (ILucas et al. 
brown dwarfs whose masses are below the deuterium-burning limit of 13 Jupiter masses 
(Mj up ), and their existence imposes severe restrictions on brown dwarf formation models. 
The models must account both for the generation of low-mass cloud fragments and also for 
their subsequent collapse. In principle, the required ma sses can be produced by mechanisms 



such as fragmenta tion of collapsing 



protostellar disks (jShen &: Wadsley 



proto stellar cores (JBoss 



200 ll ) or encounters between 



20061 ). but this still leaves the problem of providing 



the nece ssary cooling to allow grav i ty to overcome gas pressure. In this regard, recent 



work by IWhitworth fc Stamatellosl (j2006l ) suggests that masses as low as 1-4 Mj up are 
possible, the lower limit being independent of the particular star formation scenario. 
However, difficulties still arise when subsequent accretion is taken into account, since this 
can greatly increase the final mass of the object. Subsequent dynamical ejection may b e 



Reipurth fc Clarke! toll )). 



necessary in order to prevent this mass gain (see, for example 
A key constraint on these models could be provided by observational determination of 
the form of the initial mass function, and, in particular, whether there is a minimum 
mass for brown dwarf formation . The latter would provide a key test of the prediction of 



Whitworth fc Stamatellosl (120061 ). even though additional tests would be required in order 



to distinguish between the possible fragmentation mechanisms. 

With this in mind, we have begun a spectroscopic study of objects in the p Oph cloud 
core whose near-infrared colors are suggestive of low-mass brown dwarfs. This region is 
a particularly suitable site for the study of young brown dwarfs due to its youth, high 
rate of low-mass sta r formation and its relative proximity at a distance of ~ 120-130 pc 



( iLoinard et al. 



2008 



Mamajek 



2008 



Lombardi et al. 



20081 ). We have produced a list of 
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brown dwarf candidates in a 1° x 9' region of the cloud core, based on deep integratio n 



200J) 



J, H, and K s images from the stacked calibration scans of 2MASS (ICutri et al. 
together with Spitzer IRAC images at 3.6, 4.5, 5.8 and 8.0 /im. Our selection technique is 
based on the estimated effective temperature, obtained from fits of the observed spectral 
energy distributions (SEDs) to atmospheric models (Marsh et al., in preparation). In a 
pilot program of spectroscopic observations, we have observed seven candidates from this 
list using the NIRC instument on the Keck I telescope, and now present the results. In 
particular, we report the first discovery of a planetary-mass brown dwarf in p Oph. 



2. Observations and Data Reduction 

We selected seven brown dwarf candidates whose near-infrared colors suggested 
effective temperatures less than 2000 K. These objects are listed in Table [U each is labeled 
with an ID number based on the ordering of our source extractions in the 1° x 9' 2MASS 
Deep Field; the full list will be presented i n the near future by Mar sh et al. We observed 



these objects using the NIRC instrument (IMatthews fc Soiffer 



19941 ) on the 10-m Keck I 



telescope on the night of 2009 May 27 (UT). NIRC was used in grism mode (grism grl20), 
with a slit width of 3.5 pixels, to produce low-resolution (R ~ 90) spectra in the wavelength 
range 1.5-2.4 fim. Conditions were essentially clear except for light cirrus haze, and the 
seeing was ~ 0.6". 

Targets were located on the slit by first imaging them in the K-band continuum. They 
were then observed in grism mode in sets of 30 s or 60 s exposures, using a nodding pattern 
in which the telescope was dithered at intervals of 5" along the slit. Typically, 5 such sets 
were coadded. 

The coadded images were pair-wise subtracted to remove the dark current and sky 
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background, and divided by a flat-field image to remove pixel-to-pixel gain variations. 
The latter image was derived from observations of the telescope dome after subtraction of 
the dark response. It was also used in the construction of a bad pixel mask. Correction 
for the spectral response of the instrument and absorption by the telluric bands of the 
atmosphere was then obtained from observations of AO stars Oph SI and HD 161743. The 
correction was applied by multiplying the spectrum of a program star by a 10,000 K black 
body spectrum and dividing by the observed spectrum of the AO star. In doing so, it was 
necessary to dere dden the Oph SI s pectrum to correct for the known 10 magnitudes of 



visual absorption (IGagne et al. 



20041 ). For telluric correction we also utilized an observation 
of NGC 4361, a planetary nebula whose HK spectrum is dominated by the continuum of 
the hot central star. Wavelength calibration was accomplishe d using a cubic fit to a set of 



atmospheric OH lines identified in sky observations based on 



Rousselot et al. 



fl2000h . and a 



set of nebular lines from an observation of another planetary nebula, G049. 3+88.1. 

The maximum likelihood estimate of the spectrum of each program object was 
obtained using all available observations for that object simultaneously. This necessitated 
spatial co-registration of the individual subtracted nods, accomplished by adjusting the 
image offset parallel to the slit based on maximum correlation of the signal with the 
wavelength-dependent spatial point spread function (PSF). The latter was estimated 
using observations of strong calibrators. The program star spectrum was then obtained 
from a series of maximum likelihood flux estimates as a function of wavelength using all 
observations which fell within a rectangular "fitting" window of width typically 9 pixels in 
both the spatial and spectral directions, where 1 pixel corresponds to 0.15" and ~ 0.005 /zm, 
respectively; the spectral boxcar width was, however, increased to 21 pixels for #4450, the 
object with lowest S/N. The estimation procedure involved knowledge of the measurement 
noise, estimated from adjacent strips on the image, external to the strip containing the 
signal. In order to minimize the effect of spikes due to cosmic ray hits, it was necessary 
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to use trimmed averaging during noise estimation, and 3a outlier rejection during spectral 
estimation. The estimated spectra of our seven program objects are shown by the solid 
lines in Figures [T] and [2J Also included in Figure [21 for comparison, are our spectra of two 
known T dwarfs, SDSS 1254-0122 (T2) and 2MASS 1503+2525 (T5.5). 



3. Model Fitting 



The observatio ns have been fitted to a g rid o f synthetic spectra base d on photospheric 



models described by 



Allard et al. 



(12001 



2003J) and 



Hauschildt et al. 



(119991 ). and downloaded 



from the Lyon group's websitqj- The grid covers the range T c g = 100-10,000 K in 
effective temperature and log g = 2.5-6.0 cm s -2 in surface gravity; we have assumed 
solar metallicity. The temperature range is spanned by four models (COND, SETTL, 
DUSTY and NextGen), each of which covers a particular regime with respect to dust grain 
formation; the COND model is applicable to methane dwarfs (T < 1500 K). 

For each of these models we have calculated the inverse-variance weighted sum of 
squares of residuals between the model (smoothed to the resolution of the observations) and 
the observed spectrum of a given object. During this procedure it was necessary to redden 
the model spectra since objects in the p Oph cloud are seen through a substantial amount 



of extinction; for this purpose the 



Cardelli et al. 



(119891 ) reddening law was used. For each 



object, the unknowns were therefore: T e g, logg, Ay, and the model type. Maximum 
likelihood estimates of these parameters were obtained by minimizing the mean square 
residual over the wavelength range 1.5-2.4 fim (excluding 1.7-2.0 fim. to avoid the deep 
telluric absorption bands), and the results are presented in the last four columns of Table 
[TJ The corresponding model spectra are plotted as dashed lines in Figures [TJ and [2J 



] http : //phoenix . ens-lyon . f r/Grids 
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Fig. 1. — Spectra of six of the seven candidate objects. In each case, the solid line represents 
the observed spectrum and the dotted line represents the uncertainty due to uncorrelated 
measurement noise; the dashed line represents the best-fit model spectrum. A boxcar aver- 
aging window of width 9 pixels (0.045 /im) has been applied in all cases. Peak S/N values for 
the six spectra (after smoothing), in order of object ID number, are: 248, 192, 156, 126, 87, 
and 57. Note that the pronounced dips at ~ 1.85 /xm represent an artifact due to imperfect 
correction for the deep telluric absorption bands. 
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A [/J,m] 

Fig. 2.— Spectra of object #4450 and the two known T dwarfs, SDSS 1254-0122 (T2) and 
2MASS 1503+2525 (T5.5), normalized to peak fluxes of 0.075, 3.6, and 6.2 mJy, respectively. 
For clarity, the T2 and T5.5 spectra have been displaced vertically by 0.7 and 1.4 units, 
respectively. In each case, the solid line represents the observed spectrum and the dashed line 
represents the best-fit model spectrum. The dotted line represents the RMS measurement 
noise in the #4450 spectrum. A boxcar averaging window of width 21 pixels (0.1 fim) was 
used in the estimation of the #4450 spectrum, but no averaging was used for the other 
two objects. The peak SNR of the smoothed #4450 spectrum was 13.5; the peak SNRs for 
SDS 1254-0122 and 2MASS 1503+2525 were 1188 and 1848, respectively. The approximate 
extents of the water and methane absorption bands are indicated at the top. 
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One of the objects (#4450) has a spectrum which resembles a reddened version of 
SDSS 1254-0122 (T2), and the model fitting results confirm its identity as a moderately 
cool brown dwarf. The spectrum is suggestive of a low gravity object, as evidenced by 
the steeper falloff on the short wavelength side of the if -band peak with respect to that 
of the field dwarfs in Figure [2} and the displacement of the i^-band peak to a longer 
wavelength — both of these features are consistent with the much deeper H 2 absorption 
expected in a lower gravity object. 

The maximum likelihood parameter estimates based on the observed NIRC spectrum 
are: T eff = 1600+^ K, \ogg = 4.5 1^ and Ay = 24.5 with an uncertainty of ~ ±10 mag. 
The large uncertainty in Ay is due to the relatively low S/N which allows for a fairly large 
range in model temperatures. To gain more leverage for the Ay estimate and to further 
constrain the other parameters, we repeated the model fitting incorporating the IRAC 
[3.6] and [4.5] magnitudes (15.77 ±0.13 and 15.85 ± 0.10, respectively; Marsh et al., in 
preparation). We obtained: T eff = 1400 ± 100 K, \ogg = 3.0^, Ay = 15.6 ± 2.5; the 
corresponding fit to the dereddened spectrum (± IRAC fluxes) is shown in Figure [31 



4. Discussion 



The results of model fitting indicate that: 



1. The estimated temperatures for the two known T dwarfs are consistent wit h the 



previously-published spectral types (IBurgasser et al. 



2003 



Dahn et al. 



2002). They 



also indicate surface gravities characteristic of mature brown dwarfs, and zero visual 
absorption within the error bars, as expected. These two results have therefore 
provided a useful check on our o bserving and data redu ction techniques. This is 



particularly significant in view of IBurgasser et al 



(120041 ) 's finding that model fits to 
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Fig. 3. — De-reddened spectrum of object #4450 (solid line), and the best-fit model (dashed 
line) based on T e g = 1400 K and log g = 3.0. The dotted line represents the uncertainty in 
the observed spectrum and the filled circles represent the IRAC fluxes at 3.6 and 4.5 /xm. 
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some T dwarf spectra yield anomalously low surface gravities. However that bias 
occurred only for late T dwarfs (T7 onwards), and hence there is no disagreement 
with the present results. 

2. For all except object #4450, the spectroscopically-estimated temperature is 
substantially higher than the T e g < 2000 K values inferred from our previous SED 
fits. These other six objects are most likely T Tau stars and/or reddened background 
stars. One possible exception is #2974, whose error bars place it within reach of being 
a brown dwarf near the hydrogen-burning limit (its observed spectrum is consistent 
with T c ff = 2800 K fc logq = 4.0, corresp onding to a mass as low as 0.05 M Q based 



on the isochrones of 



Baraffe et al. 



(120031 )). In our forthcoming paper (Marsh et al., 
in preparation) we discuss the likely reason for the six underestimated temperatures, 
and have used the results to improve our SED fitting procedure. 

3. Object #4450 is a low-gravity brown dwarf whose estimated temperature suggests an 
early T spectral type. This is supported by the presence of weak, broad absorption 
features on the long wavelength sides of the H and K peaks (see Figure [2]), suggestive 
of methane. Definitive identification of methane will, however, need to await future 
observations at higher spectral resolution and signal to noise ratio. 



Figure H] shows the estimated location of #4450 in the log g-T e g plane, together with 
the la contour of the model fitting s olution. Also s uperp osed on the figure is a set of 



isochrones for the COND m odel, from 



1 Myr for the p Oph cloud (jLuhman fc Rieke 



Baraffe et a 



1999 



( 200 3T). If we take the estimated age o f 



Prato et al. 



2003 



Wilking et al. 



20051 ) . 



the isochrones further constrain the parameters of this object. Specifically, the intersection 
of the la contour and the 1 Myr isochrone corresponds to T eff ~ 1400 K and log g ~ 3.3; 
such a temperature reinforces our estimate of early T spectral type. The corresponding 
mass would be ~ 2 — 3 M Jup , making #4450 a candidate for the lowest-mass object outside 
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our sola r system to have been imag e d to date. Another candidate fo r this distinction is S 



Ori 70 fjZapatero Osorio et al. 



2002 



Martin fc Zapatero Osorio 



20031 ). whose parameters 



are also indicated on Figure HI We note, however, that the identification of S Ori 70 as 
a young, low-gr a vity m ember of the a Orionis cluster has been called into question by 



Burgasser et al 



([2004]), who find that the observed spectra of that object could equally well 



be interpreted in terms of an old field T dwarf of much higher mass. 

We can estimate a distance to #4450 by comparing the apparent and absolute 
magnitudes, obtained from the observations and model, respectively. If we assume an age 
of 1 Myr, then the corresponding isochrone in Figure H] indicates a range 1360-1610 K 
of T e fj values from the spectral fits; the cor responding ra n ge in Mk, obtained by spline 



interpolation of the isochronal values from 



Baraffe et al. 



(}2003|), is then 11.64-10.72. 



Comparison with the dereddened apparent K magnitude of 15.86 ± 0.23, assuming that all 
of the uncertainties add in quadrature, then leads to a ±lcr range in distance of 68-111 
pc. Although this might suggest that #4450 is closer than p Oph, based on most distance 
es timates (120-130 pc), t his range is in fact, consistent with the recent 119 ± 6 pc estimate 



of 



Lombardi et al. 



(120081 ) when account is taken of the finite thickness of the cloud along the 



line of sight (28^g pc). Our results are therefore consistent with #4450 being located near 
the front edge of the p Oph cloud. If the object is older than 1 Myr, the estimated distance 
decreases, since older objects are less bloated and therefore must be closer to maintain the 
same flux. If, for example, we assume an age of 100 Myr (the largest age within the la 
contour of Figure HI), the distance upper limit gets reduced to ~ 86 pc. 

What evidence do we have that object #4450 is a p Oph cluster member and not 
simply a foreground or background T dwarf? We can readily exclude the former possibility 
because of the presence of substantial visual extinction. Also, an old background T dwarf 
is excluded since it would have insufficient flux, based on the discussion in the previous 
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Fig. 4. — Results of model fitting for object #4450, in the log g-T e g plane. The best fit to 
the observed spectrum is indicated by the filled circle, and the l a uncertainty i s ind icated 



Baraffe et al. 



(120031 ) . For 



by the solid contour. The dotted lines represent isochrones from 
comparison, the crosses represent t he estimated parameters of S Ori 70, another arguably 



youn g object of comparable mass (jZapatero Osorio et al. 



2002 



Martin fc Zapatero Osorio 



2003|). 
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paragraph. However, since the line of sight to # 4450 intersects the southeastern portion of 



the Upper Sco association (Ide Zeeuw et al 



19991 ). we mu st consider the l atter as an another 



possible location, which w ould imply an ag e of ~ 5 Myr (IMamajek 



pc distance of Upper Sco (IMamajek 



20081 ). However, the 145 



20081 ) is almost 3a beyond the estimated distance of 



#4450 and hence membership in p Oph is much more likely. In either case, the estimated 
mass would still be < 3 M Jup . 

Our conclusions regarding the nature of this object could be tested by future 
observations at high er spectral resolution (~ 1000), with an instrument such as NIRSPEC 



(IMcLean et al. 



19981 ). particularly if the wavelength range were extende d to include th e 



J-ban d. Measurement of key flux ratios of H2O and CH4 features (see 



McLean et al. 



(120031 )) would then facilitate more accurate spectral typing, and the gravity estimate could 
be refined using the K I doublet lines near 1.25 pm. Comparison of the resulting T e g and 
logg with the COND isochrones might then lead to a better constraint on the age which, 
in the present analysis, has been assumed to be 1 Myr based on sky location in the p Oph 
cloud. 

Our detection of this object suggests that the initial mass function for the p Oph 
star-forming region extends well into the planetary-mass regime. Object #4450 is, however, 
only one of about a dozen objects within our 1° x 9' region of study whose near-infrared 
SEDs are similarly indicative of planetary-mass brown dwarfs, and we plan to obtain 
near infrared spectra of these objects in the near future. We thereby hope to accumulate 
some statistics on the low end of the brown dwarf mass function to further our goal of 
constraining the formation models. In this regard, we also look forward to the expected 
wealth of information from the Wide-field Infrared Survey Explorer (WISE) mission, 
which will provide an inventory of low-mass cold brown dwarfs in the immediate solar 
neighborhood. 
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Table 1. Observational and model- fitted parameters. 



Object ID 


Type a 


RA 


Dec 


K 


tint 

H 


A.M. C 


[K] 


log 3 d 
[cm s — 2 ] 


A v c 
[mag] 


Mod. f 


#1307 g 


bdc 


16 27 32.89 


-24 28 11.4 


14.92 


360 


1.68 


> 4800 


indct. 


36.1 


N 


#1449 


bdc 


16 27 30.36 


-24 20 52.2 


15.67 


900 


1.40 


> 5700 


indet. 


21.4 


N 


#2403 


bdc 


16 27 21.63 


-24 32 19.2 


16.42 


900 


1.70 


> 5200 


indct. 


24.5 


N 


#2438 


bdc 


16 27 09.37 


-24 32 14.9 


16.74 


900 


1.92 


> 3500 


indet. 


32.7 


D 


#2974 


bdc 


16 27 16.74 


-24 25 39.0 


16.86 


900 


2.20 




indct. 


25.9 


D 


#3117 


bdc 


16 27 17.68 


-24 25 53.5 


15.70 


900 


2.88 


3800±200 


indet. 


24.7 


D 


#4450 


bdc 


16 27 25.35 


-24 25 37.5 


17.71 


1080 


1.43 


1400+100 


3 0+ 17 


15.6 


C 


SDSS 1254-0122 


T2 


12 54 53.90 


-01 22 47.4 


13.84 


450 


1.06 


1500+100 


5.0+0.5 


0.0 


C 


2MASS 1503+2525 


T5.5 


15 03 19.61 


+25 25 19.6 


13.96 


900 


1.06 


1100+100 


5.0+0.5 


0.7 


s 


NGC 4361 h 


PN* 


12 24 30.76 


-18 47 05.4 


14.02 


450 


1.29 










Oph Sl h 


AO 


16 26 34.17 


-24 23 28.3 


6.32 


6.1 


1.74,1.80 










HD 161743 h 


AO 


17 48 57.93 


-38 07 07.5 


7.57 


6.1 


3.09 











Note. — Spectr al types and K magnitudes for the two known T dwarfs are from 



Dahn etal 



Burgasser et al 



(2003) 



Leggett et al 



(2000), and 



I 2002). The last four columns of the table represent the results of model fitting, as described in the text. 
a Designation "bdc" means that the object was a brown dwarf candidate 
b Integration time 
c Air mass of observation 

d Designation "indet." means that the quantity was indeterminate from the model fits 
e Uncertainty in Ay is ±2.5 mag for #4450, and ±1.0 mag in all other cases 
f Bcst fit model: N=NextGen, D=DUSTY, S=SETTL, C=COND 
^Corresponds to 16273288-2428116 in the 2MASS Point Source Catalog 
h Used as calibrator 



